sum([n if sum([i for i in range(1,n/2+1) if n%i==0]) for n in range(1,10000)])


all = [0] * 28123

def d(n): return sum([i for i in range(1,n/2+1) if n%i==0])

abu=[i for i in range(1,100) if d(i)>i]

for i in p60.comb(range(6965),2):
    n = abu[i[0]]+abu[i[1]]
    if n-1<28123:
        all[n-1] = 1
for i in abu:
    if 2*i-1<28123:
        all[2*i-1] = 1

sum([index+1 for index,k in enumerate(all) if k==0])

f1=1
f2=1
term = 3
while True:
    f3=f1+f2
    if len(str(f3))==1000:
        print term
        1/0
    f1=f2
    f2=f3
    term += 1


    
